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Interplay between antiferromagnetism and superconductivity is studied by using the 3-dimensional 
nearly half-filled Hubbard model with anisotropic transfer matrices t z and t±. The phase diagrams 
are calculated for varying values of the ratio r z = t z /t± using the spin fluctuation theory within 
the fluctuation-exchange approximation. The antiferromagnetic phase around the half-filled elec- 
tron density expands while the neighboring phase of the anisotropic d a ,2_ H 2-wave superconductivity 
shrinks with increasing r z . For small r z T c decreases slowly with increasing r z . For moderate values 
of r z we find the second order transition, with lowering temperature, from the d a .2_ H 2-wave supercon- 
ducting phase to a phase where incommensurate SDW coexists with d a ,2_ !/ 2-wave superconductivity. 
Resonance peaks as were discussed previously for 2D superconductors are shown to survive in the 
d x 2_ y 2-w&ve superconducting phase of 3D systems. Soft components of the incommensurate SDW 
spin fluctuation mode grow as the coexistent phase is approached. 
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I. INTRODUCTION 

Unconventional superconductivity in strongly corre- 
lated electron systems has been one of the central issues 
in the field of condensed matter physics. It is widely 
accepted that superconductivity takes place around the 
antiferromagnetic phase for many strongly correlated sys- 
tems such as high temperature superconductors [jlj, k- 
BEDT organic compounds f|||, and some of the heavy 
fermion systems ^,^|. Naively, superconductivity in 
these systems is induced by the antiferromagnetic spin 
fluctuation [||. In fact, the spin fluctuation theory 
successfully explains d I 2„ i( 2-wave superconductivity ob- 
served in various quasi two-dimensional systems [^,^|. 
Recently, attractive experimental results were presented 
concerning the relation between antiferromagnetism and 
superconductivity. 

For example, Celn3 with a cubic crystal structure ex- 
hibits an antiferromagnetic transition with decreasing 
temperature and it undergoes a superconducting tran- 
sition under sufficient pressure 25 kbar where the super- 
conducting transition temperature is T c w 0.25 K ||. In 
relation to this system, a new heavy fermion compound 
CeRhlns has been discovered recently, in which alternat- 
ing layers Celn3 and Rhln2 are stacked along the c-axis 
and superconductivity is induced by application of hy- 
drostatic pressure [^0|. It is worthwhile to note that su- 
perconducting transition temperature T c « 2.1 K of this 
compound is considerably higher than that of Celn3 jio) . 
Adding to this fact, Celrlns and CeCoIn5 with the same 
crystal structure as CeRhlns exhibit superconductivity 



at ambient pressure 1 1 1J] , so that the crystal structure 
of CeRhlns is considered to be more favorable for su- 
perconductivity than the cubic one. These facts seem 
to indicate the importance of dimensionality for the oc- 
curence of the unconventional superconductivity. 

In another recent experiment, NQR measurement un- 
der pressure on CeCu2Si2 indicated an antiferromagnetic 
instability in the superconducting phase [ ^3|Jl4| or pos- 
sible coexistence between superconductivity and antifer- 
romagnetism. It should also be mentioned that the ob- 
servation of a neutron resonance peak such as those ob- 
served in high-T c cuprates was reported for a heavy elec- 
tron system UPd2Ai3 |l5|-[l9[|. Thus it is worth while to 
study the interplay between antiferromagnetism and su- 
perconductivity with varying degree of crystal anisotropy 
or dimensionality using the spin fluctuation theory. 

The effect of dimensionality on spin fluctuation- 
mediated superconductivity was studied previously both 
with the use of phenomenological models for spin fluc- 
tuations |2(],|2lJ and from fully microscopic calculations 
based on the Hubbard model |^2| . In these studies com- 
parisons were made between two ideal cases with simple 
square and cubic lattices. In view of recent experiments it 
seems important to study the effect of crystal anisotropy 
interpolating between ideal 2D and 3D systems. 

In this paper, we discuss the phase diagram of the three 
dimensional anisotropic Hubbard model u sing the fluc- 
tuation exchange (FLEX) approximation [p3| - [2r|| . The 
crystal anisotropy or dimensionality is controlled by a 
parameter r z which is the ratio of the out-of-plane hop- 
ping integral to the in-plane one. Notice that cases of 
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r z = and r z = 1 correspond to the square lattice and 
the cubic lattice, respectively. We restrict ourselves near 
the half-filled electron density and < r z < 1. Our re- 
sults are summarized in following four features. (1) With 
increasing r z the area of the magnetic phase in the phase 
diagram tends to expand while that of the d a: 2_j / 2-wave 
superconducting phase shrinks. (2) For a moderate value 
of r z the d x 2_ y 2-wdu\e superconducting phase undergoes 
a second order phase transition into a coexistent phase 
between d x 2_ y 2-w&ve superconductivity and incommen- 
surate spin density wave (ICSDW). (3) For the same 
value of r z the spin fluctuation frequency spectra show 
a resonance peak or ridge around the antifcrromagnetic 
wave vector, similarly to the one studied previously for 
2-dimensional systems. The latter has been interpreted 
as the spin-excitonic collective mode and was assigned to 
the neutron resonance peaks observed in Y- and Bi-based 
high-T c cuprates. (4) Low frequency components of the 
spin fluctuation spectrum corresponding to the ICSDW 
ordering vector grow significantly as the phase transition 
point is approached. 

In what follows we discuss the model Hamiltonian and 
approximation procedure in Sec. II and the results of 
calculation in Sec. III. Sec. IV is devoted for general 
discussion and summary. 



II. MODEL AND FLEX APPROXIMATION 



A. Model Hamiltonian 



The model Hamiltonian we use here is; 
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where a ko . is the creation operator for a quasi-particle 
with momentum k and spin er, U is the on-site Coulomb 
energy, m„ is the number operator for a quasi-particle 
with spin a at i-site. ek is the energy dispersion of the 
quasi-particle given by 



ek = 



-2£(cos k x -\ 
-2f r z cos k 2 



cos ky) + At' cos k x cos k y 

I A+'~ U I 1„ I 



At'r z cos /c z (cos k x + cos k y ) (2) 



where t and t' are the nearest-neighbor and the next 
nearest- neighbor hopping integrals, respectively. Here- 
after, t is used as the unit of energy. Also, we intro- 
duce a parameter r z describing the anisotropy along z- 
axis where r z = corresponds to the square lattice, 
and r z = 1 to the cubic lattice. Although the value of 
the anisotropy parameter for the next nearest-neighbor 
hopping should be generally different from that for the 
nearest-neighbor, we use the same value for them since no 
essential difference will appear by definite distinction be- 
tween them. The variation of the phase diagram with this 
anisotropy is studied using the Green's function method. 



B. FLEX Approximation 

When the interaction constant U is equal to zero, the 
one-particle Green's function of the quasi-particle at tem- 
perature T is given by 



G (0) (k,i W „) 
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(3) 



where w„ = (2n + l)nT is the Fermion Matsubara fre- 
quency and /i is the chemical potential. 

For the three dimensional system with U ^ 0, various 
ordered states appear at respective transition tempera- 
ture below which features of the system is described by 
corresponding order parameters, and anomalous Green's 
functions. In order to study the interplay between mag- 
netism and superconductivity, we wish to construct the 
phase diagram of this model calculating the magnetic 
and superconducting transition points. In the FLEX ap- 
proximation the superconducting T c is calculated within 
a mean field level while the magnetic susceptibility is 
renormalized self-consistently and thus there is no Neel 
temperature Tn for 2D-systems. Here the finite value 
of r z or 3D-character makes it possible to find magnetic 
transition point Tn from the divergence of \Q ■ 

The Dyson-Gor'kov equations for the Green's functions 
and the anomalous Green's functions are given by 

G(k) = G (0) (fc) + G (0) (fc)S (1) (fc)G(fc) 

-Gt°\k)& 2 \k)F*(k) (4) 
= Gf 0) (-fc)S (1) (-fc)F f (fc) 
+G (0) (-fc)£< 2 )(-fc)G(fc) (5) 

where the abbreviation k = (k, iu> n ) is used and G(k) is 
Green's function for U ^ 0. F'(k) is anomalous Green's 
function to describe the superconducting phase. The self- 
energies T,^(k) and E^ 2 ^(fc) are given within the FLEX 
approximation as follows ]23|-]26|]: 



zW(k)=J2V eB (q)G(k-q), 
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W(k) = -J2V S u iS (q)FHk-q), 



(6) 
(7) 



V eS (Q) = U 2 [\ X s {q) + \x c (q) - \{x s (q)+X c (q)}] (8) 
W?) = U 2 [\ X S (l) \x C {q) - \{x S {q)-X C (l)}] (9) 
with 

x s (^t4IL, x<(q) = rMn^ do) 



l-UT(q) : 



i + uT(qY 



^ {c \q) = - Et G ( fc + q)G(k)±F\k + q)F*(k)] (11) 
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where q = (q, if2 n ) and f2„ = 2nirT is the Boson Mat- 

T 
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subara frequency. Also, J2 k = J2k J2n an d -^o is the 
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number of sites. x^Hl) an< I X^Hl) are spin (charge) 
susceptibility and its irreducible part, respectively. Due 
to the analytical continuation of x^ (?) from the imagi- 
nary axis to the real axis, spectrum of spin (charge) fluc- 
tuations is provided. 

For the purpose of calculating T c we linearize this 
Dyson-Gor'kov equations with respect to F^(k) or 
SW(fe) as follows: 

G{k) = Pf^l' (12) 
F\k) = \G(k)\ 2 Z (2 \-k) (13) 

where the normal self-energy £W(fc) in eq. (12) corre- 
sponds to that of the normal state. The transition tem- 
perature for superconductivity is determined as the tem- 
perature below which the linearized equation for Y,( 2 \k) 
has a non-trivial solution. The linearized equation is 
given as 

S (2) (fc) = -£[Wfc-p)|G(p)| 2 ] s c^o£ (2) W (14) 
p 

where T c is obtained by the temperature at which the 
maximum eigenvalue becomes unity. Note that the su- 
perconducting order parameter has the same symmetry 
as£( 2 )(fc). 

C. Details of Numerical Calculation 

The FLEX calculation is numerically carried out for 
each value of r z at fixed parameter values of t' = 0.2 
and U = 6. All summations involved in the above self- 
consistent equations are calculated using FFT algorism 
for the k-space with 16 x 16 x 16 meshes in the first 
Brillouin zone and Matsubara frequency sums with suf- 
ficiently large cutoff-frequency lo c = 36. Although the 
lattice size seems to be relatively small, it is considered 
that the obtained results are quantitatively modified but 
does not qualitatively change with magnification of the 
system size. When relative error of the self-energy for 
all k and u> n becomes smaller than 10 -6 , it is assumed 
that the solution is obtained for the self-consistent equa- 
tions mentioned above. The second-order magnetic phase 
transition is given by 

t/x s (q,ia, -0) = 1 (15) 

where a spin susceptibility x s (q, i^n = 0) diverges. Since 
numerical calculation can not treat any divergences, we 
always adopt following condition as the condition for the 
magnetic transition 

C/X s (q,ia i = 0) = l-e (16) 

where e = 0.002 is used throughout these calculations; 
^X s (q>0) « 1 @. It should be noted that slight mod- 
ification for the value of e > does not bring about 



qualitative changes for our results. The superconduct- 
ing transition tempertures are determined by solving the 
gap equation (14) obtained within FLEX approximation. 
The analytical continuations of x s< " c \q) to the real axis 
are carried out by using Pade approximants. 

III. CALCULATED RESULTS 

A. Phase Diagrams 

In this section, the results of FLEX calculation for the 
model mentioned above are shown. We begin by exam- 
ining variation with the anisotropy parameter r z of the 
phase diagram in temperature T vs hole-doping concen- 
tration (1 — n) space, where n is the electron density per 
site. The calculated phase diagrams are shown in the 
Fig. 1; 1(a), 1(b), 1(c), 1(d), and 1(e) are the results 
for r z = 0.1, 0.2, 0.4, 0.6, and 0.8, respectively. Each 
symbol used in these figures has the following meaning; 
closed circles show the transition to the superconduct- 
ing phase with Big-symmetry and open circles, open 
squares, open diamonds, and open triangles correspond 
to the magnetic transition with the SDW wave vector of 
Q = (7T,7r,7r), Qi = (0.88tt, tt, tt), Q 2 = (0.75tt, tt, tt), 
and Q3 = (0.887T, 0.887T, 7r), respectively. Also, closed 
squares in Fig. 1(b) and closed diamonds in Fig. 1(c) 
show the second- order magnetic transition between the 
d x 2_ y 2-wave superconducting phase and the coexistent 
phase with the SDW wave vector of Qi and Q 2 , respec- 
tively. It should be noted that all of these phase tran- 
sitions are of the second-order. Adding to these phase 
transition points, symbols of plus, times, sharp, and as- 
terisk in Fig. 1 correspond to the crossover temperatures 
at which the peak position of % s (q, 0) changes from Q to 
Qi, Qi to Q 2 , Qi to Q3, and Q3 to Q 2 , respectively. 
For example, in the case of n = 0.85 in Fig. 1(b), the 
peak position of x s (q, 0) changes from Q to Qi with the 
decrease of temperature, and then c? ;J .2_ i) 2-wave supercon- 
ductivity appears. For more detailed change in the peak 
position we need more detailed calculation using smaller 
mesh in the q-space. In the systems with nearly half- 
filled density, the peak of X s (q, 0) always appear at the 
commensurate position q = Q for moderately high tem- 
perature region as shown in Fig. 1. This fact is consistent 
with the results of previous calculations for 2D systems 
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FIG. 1. T vs 1 — n phase diagrams for (a) r z = 0.1, (b) 
r z = 0.2, (c) r z = 0.4, (d) r z = 0.6, and (e) r z = 0.8. The ab- 
breviations used here have following meaning; SDW(Q): com- 
mensurate SDW phase, ICSDW(Q n ): incommensurate SDW 
phase with wave vector Q n , SC(Bi g ): d x 2 _ y 2-wa,ve supercon- 
ducting phase, CPM(Q): paramagnetic phase with commen- 
surate peak for x 8 (Qj 0), and ICPM(Q n ): paramagnetic phase 
with incommensurate peak at q = Q n . Also, the meanings of 
all symbols is explained in text. 



An interesting result here is that with increasing r z from 
zero or 2D limit T c first decreases but slowly. Even for 
r z = 0.4, with substantial 3D character, T c is not much 
reduced from the 2D value. Explicit plots are found in 
the inset of Fig. 2. For still increasing value of r z , T c 
decreases rapidly. These results indicate that materials 
with layered structure such as cuprates and Ce-115 heavy 
fermion compounds are quite favorable for the appear- 
ance of c? a .2_„2-wave superconductivity. 
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FIG. 2. Calculated results of oj-dependences of Imx s (w) for 
respective r z -value at fixed electron density n — 0.77 where 
temperature is set just above T c except for the case of r z = 0.8 
for which T = 0.005. Inset: r z -dependence of T c at fixed 
electron densities n = 0.81 and n — 0.77. 

Another important result is the possible coexis- 
tent phase between d x 2_ y 2--wave superconductivity and 
ICSDW as seen in Fig. 1(b) and 1(c). It should be noted 
that the SDW ordering vectors Qi and Q 2 in 1(b) and 
1(c) are the same as the ones corresponding to the peak 
positions of X s (q, 0) just above the superconducting tran- 
sition temperature. Considering this new feature of the 
phase diagram it seems worth while to look at the dynam- 
ical susceptibility around the coexistent phase. It is also 
interesting to see if the spin-excitonic collective modes 
as were found in 2D systems persist for 3D systems and 
give simultaneous explanations for the resonance peaks 
observed in high-T c cuprates and in certain heavy elec- 
tron superconductor. These problems will be discussed 
in the following subsection. 



One of the important results we find from Fig. 1 is the 
gradual suppression of the d x 2 _ y 2-wave superconducting 
phase with increasing r z or increasing three dimension- 
ality in contrast to the gradual expansion of the antifer- 
romagnetic phase. For r z = 0.8 the d x 2_ y 2--wave super- 
conducting phase is very much suppressed in accord with 
the previous calculation for the simple cubic lattice p2| . 



B. Spectra of Spin Fluctuations 

In order to get some idea about the possible reason 
for the reduction of T c with increasing r z we have calcu- 
lated the q-integrated imaginary part of the dynamical 
susceptibility 
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iVo — 



iS) 



(17) 



for varying values of r z at a fixed electron density 0.77 
and just above T c . The results are shown in Fig. 2. The 
r z -dependence of T c is also shown in the inset for n = 0.77 
and n = 0.81. 

We see from this figure that with increasing r z the rel- 
atively low frequency components (but much higher than 
the order of T c ) decreases while the higher frequency side 
of the spectrum is fixed. This decrease of the intensity 
of spin fluctuation is particularly significant for r z > 0.4 
and is considered to be responsible for the reduction of 
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FIG. 3. Calculated results of (^-dependences of (upper) 
X s (q, uj + iS) and (lower) Imx s (q, u> + i(5) for q = Q2 (solid 
line) and q = Q (dotted line) at two temperatures (left) 
T = 0.0137 and (right) T = 0.0057 below T c = 0.0157 where 
parameter set (r z ,n) is fixed to (0.4, 0.83). 

We next discuss the behaviors of the predominant com- 
ponents of the dynamical susceptibility in the situations 
where the phase transition between the d x 2_ y 2--wa,\e su- 
perconducting and the coexistent phases takes place as 
seen in Fig. 1(b) and 1(c). As a representative exam- 
ple we choose the parameter set (r z ,n) = (0.4,0.83) in 
Fig. 1 (c) . We show in Fig. 3 the spin fluctuation spectra 
below T c for the antiferromagnetic component with the 
wave vector Q and for the SDW components with the 
incommensurate wave vector Q2. The upper and lower 
pannels show the calculated results of the w-dependences 
of the irreducible spin susceptibilities and the imaginary 
part of the spin susceptibilities, respectively, where these 
quantities have following relation according to equation 
(10) 



Imx s (q,w + W) 



Imx s (q, uj + 18) 



[1 - £/Rex s (q, w + i<5)] 2 + [t/Imx s (q, u + i8)} 2 ' ^ 
We first note that the intensity of Imx s (Q, uj + iS) around 



uj = 0.08 develops strongly with decreasing temperature. 
As proposed in previous papers [^8|-^0|, it is considered 
that this peak corresponds to the spin-excitonic collective 
mode within the d x 2_ y 2-wave superconducting phase tak- 
ing into account that the excitation energy for this mode 
is less than the superconducting gap 2A « 0.17 estimated 
from the peak-position of Imx s (Q, oj + 18). Noting that 
relatively large value of r z is used in Fig. 3, it is con- 
sidered that such a collective mode is not specific to 2D 
systems but appears more generally in 3D systems with 
the d ;c 2_ y 2-wave superconducting phase induced by the 
antiferromagnetic spin fluctuation. Thus, although this 
feature is discussed in relation with the 41meV resonance 
peak observed in Y- and Bi-based cuprates pi] , p2] |, it is 
expected that such a collective mode within the d x 2_ y 2- 
wave superconducting state should be seen in another un- 
conventional superconductors such as UPd2Al3 [p~5| p~0|| - 
As is seen in Fig. 3, the spectrum of the spin fluctuation 
with the wave vector Q2 also exhibits such a collective 
mode around oj = 0.02. In addition to this, a new excita- 
tion appears with decreasing temperature as a shoulder 
of Imx s (Q2, w + i£) near uj = . This behavior is consid- 
ered to be associated with the instability of the ICSDW 
mode Q2 at the second order transition point to the co- 
existent phase. Recent higher-resolution inelastic neu- 
tron scattering experiment shows that the shoulder on 
low excitation energy is observed at the incommensurate 
position below T c in the under-doped high-T c cuprates 
p3[ . Such an experimental data may be related with the 
behavior shown in Fig. 3. 
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FIG. 4. Calculated results of T-dependences of 1/TiT for 
r z — 0.1 (dotted line) and r z = 0.2 (solid line) at fixed elec- 
tron density n = 0.87 where their d x 2_ y 2-wa,ve superconduct- 
ing transition temperatures are almost same as seen in Fig. 
1. 

These low energy excitations are considered to affect 
various physical quantities, especially, the nuclear spin- 
lattice relaxation rate 1 jT\ which is expressed as 
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Imx s (q, u> + id) 



to 



(19) 



where the hyperfine coupling constant is assumed to be 
unity. It is well known that \jT\T exhibits the T 2 -law 
for the d x 2_ y 2-w&ve superconductor due to the existence 
of line nodes as in the high-T c cuprates. But, it is ex- 
pected that if this type of excitations exist as seen for 
Inr\; s (Q2, lu + iS) in Fig. 3, \jT\T should show an upturn 
below T c as shown in Fig. 4. In fact, such a behavior is 
recently reported for CeCu2Si2 fl||,|l4| where it is known 
that the unconventional superconductiviting transition 
takes place at low temperature and the magnetic phase, 
so-called A-phase, exists under the magnetic field. 



IV. DISCUSSION AND SUMMARY 

Since the second-order phase transition accompanies 
some symmetry-breaking, it is instrucitve to discuss 
phase transitions considering the symmetry for the re- 
spective phase in those phase diagrams. As well- 
known, the Hubbard model shown in cq. (1) has 
L*4; l xS'[/(2)x0xt r (l)-symmetry where D±h is tetrago- 
nal point group of the system, SU(2) rotational invari- 
ance in the spin-space, time-reversal symmetry, and 
U(l) gauge symmetry. When the antiferromagnetic or 
SDW transition takes place with the magnetic moment 
along z-axis, the SU(2) rotational invariance is broken 
and Z?4/jX9 reduces to its subgroup D^h, so that the 
symmetry of the system lowers to D 4^x11(1) where G is 
a magnetic point group which includes some group ele- 
ments for G's typical element R to couple with the time- 
reversal operator 9. Also, the d :c 2_ y 2-wave superconduct- 
ing transition (belonging to B\ g irreducible representa- 
tion) breaks U(l) gauge symmetry and 7r/2-rotational 
invariance, so that the symmetry of this phase becomes 
D 2h xSU{2)xQ. 

Of course, the first-order phase transition between 
SDW and d x 2 _ y 2 -wave superconducting phases is pos- 
sible where a group of one phase is not a subgroup of the 
other. Noting that no magnetic transition takes place in 
the d a ,2_j / 2-wave superconducting phase in Fig. 1(a), it 
seems that such first-order phase transitions take place 
between SDW and ^2^2 -wave superconducting phases 
as schematically shown in Fig. 5(a). On the other hand, 
the second-order phase transitions are expected between 
the coexistent phase and the ICSDW and d,. 2^2- wave su- 
perconducting phases because the group D 2 h of the coex- 
istent phase is a subgroup of both of the other two phases 
as seen in Fig. 5(b). It seems that Fig. 1(b) corresponds 
to such a case. In Fig. 1(c), it seems that the transition 
between the d x 2_ y 2-w&ve superconducting phase and the 
coexistent phase is of the second-order while the first- 
order phase transition is expected between the d x 2_ y 2- 
wave superconducting and the ICSDW phases, as shown 
in Fig. 5(c). 



Finally, it may be worth while to make a comment on 
the spin excitations in the coexistent phase. In view of 
the symmetry of the coexistent phase we may naturally 
expect to have a Goldstone or zero frequecncy spin wave 
mode which recovers the broken 5'?7(2)-symmetry. We 
have seen in Fig. 3 that the soft component as a shoul- 
der of the spin excitonic collective mode at the ICSDW 
wave vector Q2 grows while the main peak position is 
kept finite as the coexistent phase is approached. The 
softening component may be regarded as a precursor to 
the appearance of the Goldstone mode. Thus in the 
coexistent phase we may generally expect to have two 
types of spin excitation modes, the spin excitonic collec- 
tive mode associated with the superconducting gap and 
the spin wave mode associated with antiferromagnetism 
(including SDW). When Tn is substantially higher than 
T c , however, the appearance of the former mode itself 
seems to be dubious from consideration of the possible 
gap structure. It should be an interesting future subject 
to study the dispersions and the ranges of appearance of 
these modes. 
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FIG. 5. Schematic phase diagrams corresponding to Fig. 
1(a), 1(b), and 1(c). The left edges of these figures corre- 
spond to half filled density. Here, PM, SDW, ICSDW, and 
SC(Big) stand for the paramagnetic metal, SDW, incommen- 
surate SDW, and d x 2_ y 2-wa.ve superconducting phase, respec- 
tively. In SC(£>ig)+ICSDW phase, d x 2 _ y 2-wave superconduc- 
tivity and ICSDW coexist. Also, the symmetry of respective 
phase is written in order to discuss what kind of symmetry is 
broken. 

In summary, the phase diagram of the three- 
dimensional anisotropic Hubbard model with varying ra- 
tio r z of the transfer matrices along z-axis and in xy- 
plane is calculated by using the spin fluctuation theory 
within the FLEX approximation. The results are summa- 
rized as follows: (1) For small r z or nearly 2-dimensional 
case we have an antiferromagnetic phase around the half- 



G 



filled electron density and a d x 2_ y 2-wave superconducting 
phase neighboring to it. With increasing r z the antifer- 
romagnetic phase tends to expand while the d x 2_ y 2-wave 
superconducting phase shrinks. For relatively small val- 
ues of r z , however, this trend is weak and the supercon- 
ducting transition temperature T c decreases slowly with 
increasing r z . As r z approaches 1, T c decreases rapidly. 
This behavior of T c seems to be related with the reduc- 
tion in the q-integrated intensity of the spin fluctuation 
spectrum mainly in relatively low frequency (but much 
higher than the order of T c ) side. (2) For a moderate 
value of r z the d x 2_ y 2-wa,ve superconducting phase un- 
dergoes a second-order transition, with decreasing tem- 
perature, into a phase where incommensurate spin den- 
sity wave (ICSDW) coexists with d x 2_ y 2-wa,ve supercon- 
ductivity. (3) For the same value of r z the spin fluctu- 
ation frequency spectra show a resonance peak or ridge 
around the antiferromagnetic wave vector, similarly to 
the one studied previously for two-dimensional systems. 
These spin-excitonic collective modes are considered to 
explain the neutron resonance peaks observed in Y- and 
Bi-based high-T c cuprates and in a heavy Fermion sys- 
tem UPd2Al3. (4) The low frequency components of the 
ICSDW mode of spin fluctuations grow significantly as 
the transition point between the d x 2_ y 2-w&ve supercon- 
ducting and coexistent phases is approached. 

Recently, in order to discuss the effect of orbital de- 
generacy to superconductivity, the weak coupling theory 
is developed where the orbital fluctuation plays an im- 
portant role as well as spin fluctuation p4|j3q 1 . Since any 
fluctuations should be treated dynamically, it is desire- 
ble to construct the strong coupling theory including the 
orbital degree of freedom. 
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